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ABSTRACT 
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mittency and volatility of wind power generation, which is called Security-Constrained Wind Hydro- 
thermal Coordination (WHTC). In restructured power systems, Independent System Operators (ISOs) 
execute the Security-Constrained Unit Commitment (SCUC) program to plan a secure and economical 
hourly generation schedule for the daily/weekly-ahead market. The objective of security-constrained 
Keywords: daily hydrothermal generation scheduling is to determine an optimum schedule of generating units for 
Wind hydrothermal scheduling _ minimizing the cost of supplying energy and ancillary services with considering network security 
a Uncomment constraints. The problem formulation includes dynamic ramp-rate constraints for generation schedules 
ARIMA processes and reserve activation, and minimum up-time and down-time of conventional units. Of particular 
Stochastic programming interest in this study are considering more practical constraints and rigorous modeling of thermal and 
Load uncertainty hydro units such as prohibited operating zones and valve loading effects. Furthermore, for the hydro 

plants, multi performance curve with spillage and time delay between reservoirs are considered. To 
assess the efficiency and powerful performance of mentioned method, a typical case study based on 
modified IEEE-118 bus system is investigated and the results are compared to each other in different test 


system. 
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1. Introduction 


The electric energy securance issue has been getting noticed in the 
governments in whole of the world. Especially in recent years that the 
energy consuming has been increased due to industry deploy and 
dramatic use of electric devices instead of mechanic types. Recently, 
many of governments are decided to invest in the new resources 
particularly sustainable energy resources to providing the required 
energy. The main reasons of this plan are summarized as (a) as to 
being competed the fossil energy resources and raising its price in 
international energy market and (b) the legislation of environmental 
laws such as stated in the “Low Carbon Transition Plan”, July 2009, the 
governments are obliged to reduce the Carbon emission to 80% by 
2050 [1]. Videlicet, wind turbine generators (WTG) will supply 20% of 
the US power demand by 2030 [2] and European Union planned for 
20% power generation from renewable energy sources by 2020 [3]. 
But, the increment in penetration of wind power plants on power 
systems poses well-documented technical challenges for the existing 
power systems because of uncertainty production characteristic. 
Especially, in modern systems so-called Smart Grids, it is necessary 
to consider the unit commitment (UC) and economic dispatch 
incorporating wind farms generation. 

Therefore, uncertainty of renewable energies must be included 
as an aim function in system scheming daily and weekly means an 
economic and secure schematization which consists of the hydro 
and thermal units incorporating wind power plants. Daily Hydro- 
thermal Generation Scheduling (DHGS) is a unit commitment and 
economic dispatch among hydro and thermal units during a 
scheduling period of time for supplying the demand in minimum 
total cost. Also to obtain the realistic results, the practical con- 
straints related to thermal plants, hydroelectric system and elec- 
trical power system (i.e. satisfying power system demand) should 
be investigated. On the other hand, ISO coordinates market 
activities for satisfying the hourly load demand, limited fuel and 
other resources, environmental constraints, and transmission 
security requirements [4-6]. 

Several literatures are studied on methods to estimation of 
accurate wind speed and consequent power output of wind farms. 
These investigations have been based on such foundations as fuzzy 
logic [7], neural networks [8], and time series [9]. In this paper, the 
time series model of Auto Regressive Integrated Moving Average 
(ARIMA) has been used for modeling of wind power generation 
uncertainty [10]. 

Providing the Security-Constrained Hydrothermal Coordination 
(SCHTC) with the lowest cost is the nonlinear problem with multiple 
minima. There are some techniques such as Lagrangian Relaxation (LR) 
[4], Dynamic Programming (DP) [5], Mixed Integer Programming 
(MIP) [11], Benders Decomposition (BD) [12] and various intelligent 
techniques [13-15] to find the best solution for a goal function with 
complex and nonlinear characteristics and heavy equality and inequal- 
ity constraints. In this paper, the MIP approach is applied for solving 
the SCHTC with volatile wind power generation problem which many 
of the practical constraints for detailed modeling of the all of type 
units are considered. Ref. [16] develops a two-stage stochastic mixed- 


quadratic programming problem has been proposed a new modeling 
of the optimal unit commitment and sale bid to the day-ahead market, 
and the optimal economic dispatch of the bilateral contracts for all the 
thermal and combined cycle units. This paper has been used from real 
data of a Spanish generation company and market. 

In electricity market, the ISO has an important role in manage- 
ment of the whole network by possession of the commit and 
dispatch of all energy resources based to required load with 
meeting of security constraints and minimum cost [17-19]. Actu- 
ally, the ISO determine an optimum schedule of generation units 
with the standard market design (SMD) according to security- 
constrained unit commitment (SCUC), for minimizing the cost 
[20]. This is done while in several studies the GENCO effort to 
scheduling the generation of units to obtain the maximum profit 
[21]. In this paper the unit commitment and scheduling of system 
units (hydro, thermal and wind) is done and analyzed from ISO's 
viewpoint for having the minimum cost. 

Finally, in this paper the both thermal and hydro subsystems 
are considered in planning of ISO against of [20-23] that studied 
them separately without considering of security constraints. Also 
limitations of the ramp rate, frequent switching, frequent cycling 
of hydro units in daily operations according to mechanics pro- 
blems and costs related to start-up and wear and tear of those, and 
prohibited operating zones (POZs) in thermal units are considered 
in UC and ED [24,25]. In addition to prior mentions, the hydro and 
thermal units operating characteristics such as valve point loading 
and minimum up and down time constraints are added to the 
proposed model for more precision although may be caused to a 
higher nonlinear, non-smooth and non-convex function [24]. Also 
to conducting the simulation, the load is not constant and has 
been considered as a variable parameter. The load uncertainty is 
modeled based on the error of load prediction and needed 
scenarios for solving the problem are generated by the Monte 
Carlo Simulation (MCS) method and roulette wheel mechanism 
[26,27]. 

The remaining sections of this paper are organized as follows. 
In Section 2, the stochastic model of hybrid wind-hydrothermal 
scheduling is formulated as a 0/1 mixed-integer linear program- 
ming problem. The modeling approaches used in this paper to 
represent the load and wind power generation uncertainties are 
described in Section 3. Obtained results from case studies with 
detailed discussion are presented in Section 4, and at the end the 
conclusions are drawn in Section 5. 


2. MIP formulation of security-constrained wind 
hydrothermal coordination 


Mathematically security-constrained wind hydrothermal coor- 
dination is a decision problem with an objective to be minimized 
with respect to a series of prevailing equality and inequality 
constraints. There are many parameters, equality and inequality 
constrains for reality optimum scheduling of units in power 
system. These limitations must be considered in the goal function 
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Nomenclature 

Indices 

I thermal unit index 

J hydro unit index 

wf wind farm index 
time interval (hour) index 

s scenario index 

u uncertain parameter index 

Constants 

q conversion factor equal to 3.6 x107°(Hm? s/m? h) 
(1 Hm?=10° m?) 

O number of periods of the planning horizon (24 h) 


0g.) minimum water discharge of unit j at hour t (m?/s) 
O(j, t) maximum water discharge of unit j at hour t (m?/s) 


Tij time delay between reservoir of plant i and reservoir 
of plant j (h) 

Ai shut down cost of unit i ($) 

Aj start-up cost of unit j ($) 

b, (i) slope of block n of fuel cost curve of unit i($/MWh) 

b„(j) slope of the volume block n of the reservoir associated 
to unit j (m?/s/Hm?) 

bE) slope of the block n of the performance curve of k unit 


j (MW/m?/s) 
be,(i) slope of segment n in emission curve of unit i 


DT(i) minimum down time of unit i (h) 

ei valve loading coefficient 

€min(i) | Minimum emission generation of unit i (lbs) 
EGR emission group consists of SOz, NOx 


E(p_,() emission generation of (n—1) th upper limit in emis- 
sion curve of unit i 

EMmax(i) maximum emission of unit i (Ibs) 

EMmax(gr) maximum emission of group (Ibs) 

fi valve loading coefficient 

FGR fuel group consists of coal, gas and oil 

Fmin(i) minimum fuel consumption of unit i (Mbtu) 

Fmax(i) maximum fuel consumption of unit i (Mbtu) 

Fmin(gr) minimum fuel consumption of group (Mbtu) 

Fmax(gr) maximum fuel consumption of group (Mbtu) 

F(p}_,(i)) cost of generation of (n—1)th upper limit in fuel cost 
of unit i 

F(j,t,s) forecasted natural water inflow of the reservoir asso- 
ciated to plant j in period t (Hm?/h) 

KG) Cost of the Ath discrete interval of the start-up cost of 
unit i ($/h) 


I(i) initial status of unit i (0/1) 

L number of variable head 

M number of prohibited operation zones 

NS number of scenario after scenario reduction 
p(s) probability of scenario s 


p ™°™(s) normalized probability of selected scenario s 
MSR(i) maximum sustained ramp rate (MW/min) 


MU maximum number of units that can be on at the 
same time 

NB number of bilateral contract 

NES number of discrete intervals of the start-up emission 
function 

NE number of discrete intervals of emission function 

NL number of blocks of the piecewise linearization of the 


variable cost function 


AC) power capacity of bilateral contract m at hour t (MW) 
Pmin(i) minimum power output of unit i (MW) 
Pmax(i) maximum power output of unit i (MW) 


Pp, minimum power output of plant j for performance 
curve n(MW) 
DU) capacity of plant j (MW) 


psi) lower limit of nth prohibited zone of unit i (MW) 
pł (Ò upper limit of (n—1) th prohibited zone of unit i (MW) 


Q) minimum water discharge of hydro plant j if is on (m?/ 
s) 
QnG) maximum water discharge of block n of plant j (Hm?) 


QSC(i) quick start capacity of unit i (MW) 
i 


RDL,(i) ramp down limit for block n (MW) 

RUL,(i) ramp up limit for block n (MW) 

s°(i) time periods of unit i has been shut-down at the 
beginning of the planning horizon (h) 

S(j) maximum spillage of unit j (m°) 


Smax(i) maximum hour unit i can be off (h) 


SUE(i) | emission generated by unit i when started up (lbs) 
SDE(i) | emission generated by unit i when shut down (lbs) 
SDF(i) fuel consumption of unit i when shut down (Mbtu) 
SD(i shut-down ramp rate limit of unit i (MW/h) 


) 

) start-up ramp rate limit of unit i (MW/h) 
UT(i) minimum up time of unit i (h) 

i) 


Ui time periods of unit i have been on-line at the 
beginning of the planning horizon (h) 

volj) minimum content of the reservoir associated to plant j 
(Hm?) 

v°(j) reservoir content at the beginning of the study time 
(Hm?) 

vecj) reservoir content at the end of the study time (Hm?) 

Vn(j) maximum content of the reservoir j associated to nth 


variable head (Hm?) 
Variables 


6n(i,t) generation of block n of fuel cost curve of unit i at 
hour t 

y(i, t) dummy variable (h) 

y n(i,t) generation of block n of unit i at hour t of valve point 
loadings curve 

Wn volume block n for the reservoir of hydro plant j at 


hour t (MW) 

B(i,t) start-up cost of unit i at hour t ($) 

b,(i) slope of power block n of fuel cost curve of unit i 
($/MWh) 

biG) slope of the block n of the performance curve | of 
hydro plant j (MW/m?/s) 

c(i,t) valve point loadings cost of unit i at hour t ($) 


F(i,t) fuel cost of unit i at hour t ($) 

Nit) non-spinning reserve of a unit i at hour t for bid on 
spot market when unit is off (MW) 

N4(j,t) non-spinning reserve of a unit j at hour t for bid on 
spot market when unit is off (MW) 

N,(it) non-spinning reserve of a unit i at hour t for bid on 
spot market when unit is on (MW) 

N,G,t) non-spinning reserve of a unit j at hour t for bid on 
spot market when unit is on (MW) 

plit) real power generation of unit i at hour t (MW) 

Pmin(1. t, S), Pmax(i t S) lower and upper limit of real power 
generation of unit i at hour t (MW) 

pt) real power generation of unit j at hour t (MW) 
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pG,t,s),pG,t,s) lower and upper limit of real power generation 
of unit j at hour t (MW) 
p(t) power for bid on spot market at hour t (MW) 


Qt) water discharge of unit j at hour t (m?/s) 

qn(jt) water discharge of block n of unit j at hour t(m?/s) 

R(i,t) spinning reserve of a unit i at hour t for bid on spot 
market (MW) 

RG,t) spinning reserve of a unit j at hour t for bid on spot 


market (MW) 
RDL(p(i,t)) ramping down limit of unit i at hour t (MW) 
RUL(p(i,t)) ramping up limit of unit i at hour t (MW) 


(i,t) time periods that unit i has been shut-down at hour 
t (h) 

s(j,t) spillage of the reservoir associated to unit j at hour t 
(m?/s) 

v(j,t) water content of the reservoir associated to plant unit 


j at hour t (Hm?) 
Binary variables 


f,(i,t,s) 1if block n of fuel cost curve of unit i at hour t selected 
£,G,t,s) 1 if variable head n+1 of unit j at hour t selected 
Xn. t,s) 1 if power output of unit i at hour t has exceeded block 


n 

h,(j,ts) 1 if the water discharge of unit j at hour t has exceeded 
block n 

Kits) 1 if unit iis on-line at hour t 

Ij,t,s) 1 if hydro plant j is on-line at hour t 


Igi,t,s) 1 if unit i at hour t provide non-spinning reserve when 
unit is off. 

Id,(i,t,s) 1 if block n of ramping down limit curve of unit i at 
hour t selected 

Tu,(i,t,s) 1 if block n of ramping up limit curve of unit i at hour t 
selected 

w‘(i,t,s) 1 if unit iis started-up at the beginning of hour t and it 
has been offline for 2 hours 

y(it,s) 1 if unit i is started-up at the beginning of hour t 

y(j,ts) 1 if unit j is started-up at the beginning of hour t 

2(it,s) 1 if unit i is shut-down at the beginning of hour t 

z(j,sts) 1 if unit j is shut-down at the beginning of hour t 

Sets 

G set of indices of the group units 

I set of thermal units 

J set of hydro units 

WF wind farms 

N set of indices of the blocks of the piecewise lineariza- 
tion of the unit performance curve. 

T set of indices of the periods of the market time 
horizon 

S scenario numbers 

SP stochastic parameters (wf indicates wind farms) 

i,p,q,r sets of ARIMA family model 

A set of the discrete intervals of the start-up cost 
function for thermal units 

Qj set of upstream reservoirs of plant j. 

N time number of wind power generation time series. 


as security constraints included the transmission flow and bus 
voltage constraints or the start-up limitations caused extra invest- 
ment. In this section, the objective function and its different parts 
of simulation is explained clearly. 


2.1. Objective function 


In this work, the suggested objective function to determine the 
optimal planning of coordination of wind, hydro and thermal units 
over scenarios caused minimum cost is described as (1). It is 
necessary to remind that the load in this paper is varied con- 
stantly. 


Cost! = ¥ p"°™(s). cos t(s) 
S 


cos t(9= 5 [zra t, s) + AjZ(i, t,s) + Bi, t,s) + C(i,t,9] + SAG. no} 
jeJ 


teT (iel 
(1) 


The last term of variable cost function of each scenario 
represents the corresponding cost of hydro units regarding their 
start-up cost over the given period [28]. While the first term 
shows thermal operating cost including fuel, shutdown, startup 
costs and valve point loadings cost. The operating costs of wind 
units are assumed to be zero. The list of symbols is presented in 
the Nomenclature section. 

When minimizing the total cost in power systems, the total 
generation of hydro, wind and thermal plants should be equal to 
the total system demand plus the transmission network loss. Due 
to simplicity, the network loss is not considered in this paper. This 
equation is given by 


Epit, + YpG.ts9+ X pwf,t.s)=Ppe yte T, VseS (2) 
iel jeJ WF 


wf e 


2.2. Cascaded-hydro units’ model 


According to previous sections, the hydro power plants have 
the start-up limitation. Because of the unnecessary commitments, 
loss of water during start-up period, wear and tear of the windings 
and mechanical equipment, the start-up limitation must be 
considered in planning. In this regard, the relation among gener- 
ated power, water discharge, and multi performance curves of 
hydro power plant is shown in Fig. 1. The proposed model takes 
into account head dependent reservoirs with MIP formulations as 
well as hydro plants connected both in parallel and series shown 
in Fig. 2. The number of the heads for power plants is considered 
to be L in the following formulations represented in Fig. 1. 

Performance curves based on the water volume can be given by 
(3)-(6) show that the volume of each hydro power plant must be 
bigger than the minimum content of that hydro plant. Eqs. (4) and 
(5) determine the right head proportional to volume. Combination 


Generation (MW ) 
A 
BEG) Ma) SV Ut) v0) 
Pies b50) 
20) 
j bU) qt) 
P,Q) — a) 
nt) 
aG) 
vG) SvG) SvG) 
Head 2 
P,G) vod) S¥U.1) $040) 
Head 1 
P,Q) 7 
— + Water disch arge (m° / s) 
QG) QG) 


Fig. 1. Three-dimensional piecewise linear non-concave unit performance curve 
for hydro plant j at hour t. 
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Fig. 2. Hydraulic topology of the river basin. 


of 0-1 for these binary variables is prevented by Eq. (6). Bn (i,t,5) 
equals to 1 when (n+1)th head is used or water volume of the 
reservoir is greater than V,(j). In other words, these equations 
choose the right curve for head according to the content level. 


v(j,t,S)> Vo) VjeJ,VvteT,VseS (3) 
vG, t, S) < Vi@)P,_10.65) 


L 
+ È Ynn, t,S)—-Pn10.t.9] YjeJ,YyteT,YsesS (4 


n=2 
v, t, S) > V1) h,_10. t, S) 


L 
+ x, Vn-20)Bn20,0-Br10.69] VieJ,VteT, VseS (5) 


n= 


Pd t.s) > P20, t,S) > ... > Pi10, t, S) VjeJ,VteT,VseS (6) 


2.2.1. Piecewise linear formulation of the variable head multi 
performance curves 

In most papers the effect of the variation of the head has been 
neglected and head is assumed fix. The proposed MILP model take 
into consider head dependent reservoirs with a piecewise linear 
approximation. After determination of each unit in this subsection, 
the linear relationship between the generated powers, the dis- 
charged water and variable head for corresponding performance 
curve can be presented as follows: 


PU.ES)—P,DIG.t3)— X dnl t,s)bn@) 
ne 


k=1 L=1 
—PU) ley 2 Pu. fs) + Bits) <0 


n=K 


VjeJ,VteT, VseS,1<k<L (7) 
PUL SI—-PKDIGLS— E anG, t s)bh@) 
ne 


+P() 


n=k 


k=1 =A 
(k-1)— X paG, t, s) + Pabst >0 
eat 
VjeJ,WteT, VseS,1<k<L (8) 


where p(j,t,s) is the generated power by the hydro unit j at hour t, 
Pj) is the minimum power generation of the head k which is 
determined by /7,(j,t,s). Also, P(j) is the capacity of hydro unit j, and 
qn(j,t,s) is the water discharge of the block n. Finally, b,'(j) is the 
slope of the block n of the variable head k of hydro unit j. 


2.2.2. Water discharge limits 

The water discharge limit and initial and final volume and 
water balance equations same as the equations presented in [22]; 
however, in the proposed WHTC model the spillage effect is also 
considered which is inspired by [29]. 


2.3. Thermal units’ model 


In this subsection deal with the linearization of all nonlinear 
equations related to thermal units to represent the linear model 
for thermal units. 


2.3.1. Linear fuel cost function considering POZ 

The thermal units have some prohibited operating zone due to 
steam valve operation or vibration in its shaft bearing and some 
faults in the machines or their accessories such as pumps or 
boilers, etc. It has been shown that the input-output curve of 
thermal units has some valve points. These valve points generate 
some prohibited zones, so the input-output curve of thermal unit 
will be discontinuous. The quadratic production cost function can 
be accurately approximated by a set of piecewise blocks. So, the 
linear formulation of power generation for the ith thermal unit 
considering M POZs is modeled as follows: 

M+1 

F(i,t,s)= 


2 PG t, Fp- Ù) 


+bp(i)dn(i,t,s)] Viel, VteT,VseS (9) 


M+1 
Pii,t,s)= $ [PaP t, s) 
ever Viel, VteT,VseS (10) 
6n(i,t,s)>O0 n=1,2,....M+1,Viel, VteT,VseS (11) 
On(i, t. 8) < [PED- p! Dait, n=1,2,....M+1, 
Viel, VteT,VseS (12) 


M+1 
E B,(i,t.s)=Mi,t,s) Viel, VteT,VseS (13) 
n=1 


Pri, t,s) € {0, 1} 


where, po) = Pmin Ò) and Piri (0) = Pmax) 


n=1,2,...M+1,Viel, VteT, VseS (14) 


2.3.2. Cost of valve point loadings 

In practical, the generators have the multiple valves in steam 
turbines. The opening and closing of these valves caused to 
maintain the active power balance. However it adds the ripples 
in the cost function [30]. To calculate and consider the valve-point 
effects, sinusoidal functions are added to the quadratic cost 
functions as follows [31-33]: 


F,(P;) = ai + b;Pi + ¢;P? + abs(e; sin (fj(Pimin —P;))) Viel (15) 


where e; and f; are the coefficients of generator ith considering 
valve point loading effect. 

But the added sinus term to cost function cause to having a 
nonlinear and non-convex function that can be obstacle in MIP 
model. In this paper a linear formulation is proposed to overcome 
this problem as shown in Fig. 3. 


$ 2f je; ki n r 
C(i, t; S) = ax v2 È [anli t, S)—Wan+4(i, t,s)] 
n= 


ki 
+2-V2) X [Wana t, 9)-Wansald cs) Viel, VteT, VseS 
n=0 


(16) 
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> Generation (MW) 


Panli) Pini) +0/4f; Prini) +20/4f; Piali) +37/4f; Panli) 


Fig. 3. Linear approximation absolute sinus function valve point loading. 


where, Yh (i,t,s) is power generated by nth block in sth scenario. 


ki F n 
DUES) = Prin D9) + X |Van 0 tS) + Wangali ts) 


n=0 

t+Wan3(it.8) + Wanait.s)| Viel, VteT,vseS (17) 
= (i,t,s)< (i,t,s) < =i, t,s) Viel, VteT,VseS (18) 
ape! bar} <i abs ~ Af; ats > > 
Ta : T : : 
T; ni 65) < Yp t, $) < JFXn li t,S) Viel, VteT, VseS,n=2,3,...,X; 

1 1 

(19) 

Ants) €{0,1} Viel, VteT, VseS,n=1,2,...,x; (20) 


where, 
ki = floor[fi(Pmax()—Pmin(1)/7)] and Xi = floor[4f (Pmax()—Pmin(D/2)] 


Constraint (17) states that power output of unit i at hour t is the 
sum of minimum power output if unit is on, plus the power 
generated in each block. Constraint (18) limits the power gener- 
ated in first block, this power should be greater than zero and 
smaller than or equal to 7/4fi that is “power length” of each block. 
In this constraint, I(i,ts) is used ensuring that power generated of 
first block equal to zero, if unit i is off at hour t. To limit the 
generated power in each block y,,(i,t,s) is introduced in con- 
straints (18)-(20). This binary variable equal one if the generated 
power of unit i at hour t has exceeded block n. 


2.3.3. Dynamic ramping up/down limit 

Variations of load have been responded by generators by 
decreasing and increasing of generating power. But the rate of 
the power output changing is limited by two constraints. The 
increasing and decreasing constraint is defined by up and down 
ramp rate that are specified for each generator. In this work a 
dynamic ramp rate is used, as the ramp rate is a function of the 
generating power. Also in this method the M POZ are considered. 
The formulation related to dynamic ramp rate is given as 


M+1 
RUL(p(i, t,s))= > RUL,AOP,, (i, ts) Wiel,VteT,VseS (21) 
n=1 


M+1 
RDL(pii,t,s)) = X, RDLAWP,(i,t,s) Viel, VteT,VseS (22) 
n=1 


Constraints (21) and (22) indicate dynamic ramp limit with M 
POZs. The binary variables are used to indicate which operating 
zone has been selected. Another formulation for dynamic ramp 
rate is proposed in [34]. 


2.3.4, Stair-wise formulation of the time-dependent startup cost 
function 

Starting up generating unit costs money in addition to the 
running cost considered in economic dispatch. This cost is mod- 
eled as a nonlinear (exponential) function of the number of hours 
a unit has been off. In [35,36] two different methods have been 


suggested for time dependent linearization of this function. In this 
paper is used of formulation proposed in [35]. 


2.3.5. Minimum up-time (MUT) and minimum down-time (MDT) 
Another major source of complication in system planning is the 

nonlinear minimum up/down time constraints which has to be 

checked. In this study a linear formulation for these constraints is 

represented as follows: 

B(i) 


£ [0 -IG t,s)]=0 


t= 1 


viel (23) 


Tı 
È Ki, t,s)> UTG)y(i,t,s) Tı =t+ UTÒ-1, 


tT=t 

Vt=B(i)+1,...,90—UT(i) + 1 (24) 
0 

È Mi, t,s—y(i,t,]>0 Viel, Vt=O-UT;+2,...,0 (25) 
t=t 
c(i) 

È IG,t,s)=0 Viel (26) 


tsi 


T2 
È (1-1, 7, s)] > DT) Zit, Tz =t+ DT(Ò-—1, 


tT=t 


Vt=C(i)+1,...,90—DT(i) +1 (27) 


$ [1—(i,t,s)—z(i,t,s]>0 Viel, ¥t= O-DT(i)+ 2,...,0 (28) 
C=St 

Once the operation of a unit begins, the operational status should 
be maintained as least for a certain time (MUT). Eq. (23) deter- 
mines the initial status of units so unit will be on if B(i) is less than 
UT (i). The MUT limit for consecutive periods and for the last hours 
is satisfied by (24) and (25) respectively. Similarly, if the operation 
of a unit is terminated, the un-operational status should be 
maintained at least for a certain time (MDT) [36]. Constraints 
(26)-(28) enforce the MDT limit. 

The detailed formulation of other hydro units constraints such 
as initial and final volume [22], water balance [29] and operating 
services [37] considered in the problem, are represented in the 
mentioned references. 


2.4, Security constraints 


The security constraints required in security-constrained wind 
hydrothermal coordination problem, are obtained based on DC 
power flow model due to it has more accuracy than the linear 
power flow model. In this model the physical flow in a grid is 
controlled by Kirchoff's current law (KCL) and Kirchoff's voltage 
law (KVL), while in linear power flow model only the KCL is 
considered. In DC power flow model, the transmission constraints 
are considered as linear constraints. So the practical constraints 
related to thermal plants, hydroelectric system are expressed 
by piecewise linear approximation, therefore the security-con- 
strained WHTC is a Mixed Integer Linear Programming (MILP) 
problem. The related formulation to the DC power flow is given as 


NG, Lp 
E Pe-PP = X Pe Yb,vyteT (29) 
i=1 [=I 
1 
Fe = (Os—Sy) vb,yteT (30) 
l 


Transmission flow limits in the base case: 


FR < Fy < Fp% Vl, YteT (31) 
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3. Wind power generation and demand forecast uncertainty 


One of the most important issues considered recently by 
researchers of restructure power systems widely, is the uncertainty 
in operation of power systems. Uncertainties such as uncertainty 
due to system load prediction error or uncertainty due to unex- 
pected departure of network equipment (i.e. production units and 
lines). One another these uncertainties related to the error in 
prediction of the accurate value of wind farms output, especially 
which penetration of this energy resource type in power grids are 
increasing. 


3.1. Wind power uncertainty modeling 


The wind power is a non-predictable complete accurately. This 
future is rooted from the wind nature. Therefore the primary 
problem for incorporating wind power plants in power system is 
the uncertainty of them that decrease the reliability of grid. 
Several investigations have looked at the prediction of wind speed 
for use in determining the available wind power [7,8]. 

In literatures that study on the wind farms uncertainty, the 
wind power plant is modeled on two method, wind speed-based 
approaches [38-42], and wind power-based approach [10]. In the 
wind speed-based approaches, the wind speed must be measured 
accurately and the generating power be calculated with model of 
wind-turbine that usually is unavailable. While in power-based 
approach does not require the model of wind farm due to the 
output of them is measured and recorded directly, so the whole of 
wind power data with its details is available. Also in this method 
the error of prediction is smaller due to the amount of generated 
power by a wind-turbine is proportional by the cube of wind 
speed, so an small error (€) in the wind prediction produce a more 
big error (e°) in expected power output than when the wind 
power is measured directly. 

Although there are other methods for consideration of volatility 
of the wind power generation like as Monte Carlo simulation, and 
Latin hypercube sampling method [43], in this paper, inspired by 
{44], an Auto-Regressive Integrated Moving Average (ARIMA) 
process is used to model wind power generation uncertainty 
directly as a wind power-based stochastic approach. In the 
following, the ARIMA process is described clearly. 


3.1.1. ARIMA model of wind power 
An ARIMA (p,d,q) model of the non-stationary random process 
W(t) is depicted as [10]: 


(- 3 oo’) (a—Dy' W(t) = 99 + (: £ aw) n(t) (32) 
= i 


(=l 

where {gj} is the Auto-Regressive (AR) coefficients; {0;} is the 
moving average (MA) coefficients; n(t) is a white Gaussian process 
with zero mean and variance o2; the parameter Qo refers to the 
deterministic trend term when r>0. 

The coefficients of the ARIMA model can be estimated by 
different approaches, e.g. the Yule-Walker estimator, the least 
square estimator, and the maximum likelihood estimator [45]. For 
specification of the ARIMA's family there are two model selection 
criteria in various stochastic analyses that the Akaike's Information 
Criterion (AIC) [45] and Bayesian Information Criterion (BIC) is 
more common than other types. The ARIMA(0,1,1) model that is 
appropriate for the wind power time series is described as Eq. (33). 


(1—D)Wo(t) = (1-4 D)a(t) 


WP(t)=W32(t) for t=1,...,N (33) 


In this paper, the ARIMA(0,1,1) are used in wind-hydrothermal 
coordination for generating of different wind power generation 
scenarios and modeling of random trend of the wind farm. 

The relation between the input wind power and the output 
electric power relies on several factors, such as the efficiencies of 
generator, wind rotor, gearbox, and inverter, depending on what 
kind of turbine under investigated. For a generic wind turbine, 
some researchers [46,47] used a simplified model to characterize 
the relation between the wind power generation (WPG) and wind 
speed [48]: 


w=0 for v<v; and v>Vo (34) 

w=w r for Vi<V<Vo (35) 
Vr—Vi 

w=w, for Vr<V<Vo (36) 


Eqs. (34)-(36) show the power output of wind power plants is 
limited by wind speed, while the ARIMA model cannot consider 
these limitations. Thus these margins has been considered for 
output wind farms and added to the standard ARIMA model to 
include the upper and lower bounds of the WPG as follows: 


Wimax Wolt) > Wmax 
Wolt) = ¢ Wott) Wmin < Wolt) < W max (37) 
Wmin Wolt) < W min 


where Wmax and Wmin denote the upper and lower bounds of the 
square root of the wind power output, respectively. The block 
diagram of the limited ARIMA(0,1,1) is shown in Fig. 4. 


3.2. Load uncertainty modeling 


In this paper, the load uncertainty is considered as the error of 
load prediction. Thus, the probability distribution function of load 
prediction error must be derived based on the records of earlier 
loads. It should be mentioned; here it is assumed that the load 
forecasting process as one of the tasks of the ISO, before the 
implementation market has been done. Another noticeable point 
is that, considering each of the loads connected to the network as 
an independent variable will be caused to increasing of the 
number of random variables and the problem will be more 
complex. To solve this problem, the network load is considered 
as a random variable. It is natural that for getting the value of each 
of the loads connected to the network, it is sufficient to be divided 
based on the ratio of each bus load to the total load in the base 
state. Fig. 5 is an example of the probability distribution function 
of the prediction error of network load and with its discretization. 

As is seen in Fig. 5, seven different load intervals centered mean 
error of zero (base state) is considered, so that distance of between 
the different levels of load forecast error equal to the standard 
deviation (SD) of load prediction error. 

The important issue is how to model the random load level as a 
random variable. To implement this is used from the roulette 
wheel mechanism [49,50]. For this aim, first is applied the 
possibility of different levels of load prediction error based on 
the unit load so that the total probability is equal to one. Then, 
the distance between range [0,1] is covered based on the normal 
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Fig. 4. Block diagrams of ARIMA(0,1,1)-type wind power time series models [21]. 
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Fig. 5. Typical discretization of the probability distribution of the load error 
[26,27]. 
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Fig. 6. Roulette wheel mechanism for the normalized probabilities of the load 
forecast levels. 


probability of forecast error levels of load. Naturally, whatever the 
probability of prediction error levels of load is more; it will occupy 
more space of the roulette wheel. To clarify this mechanism, the 
roulette wheel in Fig. 6 is illustrated. 

After forming of the roulette wheel, a random number in the 
range [0,1] is generated. The number produced will be in specified 
ranges on the roulette wheel are associated with different levels of 
load forecast error, this means that the load level related to this 
interval are selected for the desired scenario. Thus, for generating 
of each new scenario, the roulette wheel mechanism will be used 
to model the behavior of random load. 

A higher number of scenarios results in a better modeling of 
uncertainties but with the cost of higher computation burden. 
Accordingly, after the production scenarios, it may be a very low 
probability of some scenarios and will only cause to increasing of 
time and volume calculations. Therefore, it is necessary to elim- 
inate the Low-value scenarios from the set of produced scenarios 
and decrease the number of scenarios. In this paper, the strategy 
adopted for this purpose is the removing scenarios with less than 
probability and the similar scenarios [51]. This reduction process 
results a model of uncertainties with more less accuracy and 
negligible error with the cost of higher computation burden. 


3.3. Scenario aggregation 


The idea of stochastic security-constrained WHTC is to con- 
struct or sample possible options for uncertain circumstances, 
solve the deterministic security-constrained wind hydrothermal 
coordination problem for the possible options, and select a good 
combination of the outcomes to represent the stochastic solution. 
Two methods are usually considered for scenario aggregation (the 
combination of the scenario outcomes) of the stochastic security- 
constrained wind hydrothermal coordination [52]: (1) using mini- 
mum value of the random variables over the scenarios and 
(2) applying weighted-average (expected value) of variables over 
the all scenarios. Using minimum value of daily operating cost may 


lead to pseudo-scheduling with optimistic operating costs in the 
planning horizon. Because it may be, minimum daily operating 
cost be related to the scenario with minimum occurrence prob- 
ability and has a significant difference ratio to other scenarios 
determined daily operating cost of problem. So, can be inferred 
that using of first method is a very conservative method and not 
economical. While in second method videlicet weighted-average 
or expected value, scenarios with more probability will have a 
more contribution on specification of final value of daily operating 
cost. It seems which second method is more reasonable of 
economic view point. For this reason, in this paper the second 
method is used for aggregation of different scenarios result paper 
to determine total daily operating cost of problem. 

So, in the proposed stochastic optimization framework, the 
later method (expected value) is considered. In this way, the 
solutions obtained from different scenarios are aggregated based 
on the probability laws to yield a single solution, describing the 
most probable outcome of the power system based on the 
evaluated scenarios, considered as the result of the proposed 
stochastic security-constrained wind hydrothermal coordination 
framework. 

As stated, the MCS method is implemented to simulate random 
characteristics of power systems load and then the scenario 
aggregation technology is used to convert the stochastic variables 
of the stochastic security-constrained wind hydrothermal coordi- 
nation problem into deterministic ones. A major advantage of 
scenario aggregation technique is that not only individual scenar- 
ios become simple to interpret but also the underlying problem 
structure is preserved. After running the proposed security- 
constrained wind hydrothermal coordination scheme for the 
accepted scenarios resulted from the scenario reduction, the 
results are aggregated according to the probability of scenarios 
to get the expected results of the formulation of hybrid wind- 
hydrothermal scheduling considering uncertainties. 

The aggregation is done for the scenario dependent decision 
variables I(i,t,s), I(j,t,s), F(i,t,s), p(its), pG,ts),R(@ts), RG,t,s) of the 
optimization problem. The aggregation is done as 


f= X pef, vses G8) 
ssl 

where, f is the variable that is aggregated and f, is the variable 

value at scenario s. It is noted that the objective function of the 

proposed formulation for security-constrained WHTC problem in 

Eq. (1) is also an aggregation of the objective function values of the 

scenarios. 

In the present paper, two main categories of power system 
uncertainties (load and wind power uncertainties) are considered 
in security-constrained wind hydrothermal coordination program 
that execute by ISO to plan a secure and economical hybrid wind- 
hydrothermal generation scheduling. 

When the number of variation and constraint increased by 
increasing the size of system, the solution time will increased and 
sometimes solver cannot solve the problem, but there are several 
method to overcome this obstacle. Benders decomposition approach 
[53] is useful method for this situation, also parallel processing may be 
used for solving this problem, for solving method with non linear 
constraint, linearization maybe handsome. 


4. Case studies 


A modified IEEE 118-bus test system is used to test the proposed 
algorithm for the security constrained day-ahead wind-hydrothermal 
power scheduling [54]. In addition to hydro and thermal units, two 
wind farms are assumed which are supervised financially under ISO 
operating scheme. Random behavior of wind power generation is 
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Table 1 
Results of the scenarios of the stochastic security-constrained WHTC framework. 


No. Scenario no. Normalized Total wind Daily operating cost ($) 
probability power (MW) 


1 1 0.788 1572.175 755704 
2 26 0.030 1238.43 839811 
3 28 0.013 1639.317 802221 
4 39 0.010 1557.42 742289 
5 40 0.008 1393.175 861473 
6 45 0.009 1466.166 761167 
7 52 0.010 1257.439 865732 
8 56 0.015 1497.414 867398 
9 78 0.005 1305.838 801530 
10 79 0.033 1551.793 731831 
11 82 0.013 1250.692 854589 
12 87 0.003 1565.795 770844 
13 97 0.006 1445.347 775077 
14 101 0.009 1353.77 850058 
15 105 0.024 1504.765 777797 
16 133 0.006 1504.793 789590 
17 145 0.005 1408.149 804971 
18 158 0.001 1640.804 826181 
19 164 0.006 1542.076 765741 
20 194 0.006 1232.111 895407 


modeled by ARIMA(0,1,1) time series model as described in Section 3.1. 
Two wind farms are separately modeled to simulate correlation 
concerns of the stochastic wind power behavior in different areas of 
the network. Auto and cross correlation coefficients related to two 
applied wind farms are considered as Appendix I. As mentioned in 
Section 3, the proposed MIP model of security-constrained wind 
hydrothermal coordination includes both uncertainties of wind power 
and load. The practical constraints of thermal and hydro generation 
units that detailed their modeling are mentioned in Section 2 and 
additional system-wide constraints such as fuel constraints and 
emission limits [36,38,55] and spinning and operating reserve require- 
ments [56] are considered in the stochastic optimization framework. 

The proposed model is implemented on a Pentium IV, 3 GHz 
with 1GB RAM using MILP solver CPLEX 9.0 in the GAMS environ- 
ment [57]. 

With the use of ARIMA method for modeling the stochastic 
behavior of wind power generation and roulette wheel mechanism 
to model uncertainty of load, 200 scenarios, including daily wind 
power time series and daily load profiles, are generated. It imposes a 
high computational burden to solve the security-constrained WHTC 
problem for all of these scenarios. So, the set of generated scenarios 
(200 wind power generation scenarios/ daily load profile) is reduced 
using the scenario reduction technique. The generated similar scenar- 
ios and scenarios with probability lower than 0.003 are discarded. 
Number of remaining scenarios after scenario reduction is equal to 20, 
which results in 200/20= 10 filtering ratio. So, the scenario reduction 
technique significantly reduces the computation burden of the pro- 
posed stochastic security constrained WHTC framework. At the same 
time, the most probable and dissimilar scenarios are retained while 
maintaining a good approximation of the uncertain behavior of these 
uncertainty resources. For the remaining set of scenarios, the proposed 
stochastic MILP model of security-constrained wind hydrothermal 
coordination is run. Selected scenarios, their normalized probability 
and total daily power generation of wind farms (generated power by 
WF1 and WF2) are presented in Table 1. 


Case 1. Security constrained wind-hydrothermal coordination 


In state that security constraints are considered, the minimum 
and maximum daily operating cost are related to scenario 10 and 
scenario 20 with 731831$ and 895407.38 $ respectively. The 
commitment schedules for these scenarios are shown in 


Table 2 


Scenario 10: WHTC with security. 


Hydrol 
Hydro2 
Hydro3 
Hydro4 
Hydro5 
Hydro6 
Hydro7 
Hydro8& 


Table 3 


Daily operating cost:731830.98 $ 
Hours (0-24) 


1000000000000000000000000 
1111111111111111111111111 
1111111111111111111111110 

1000000000000000000000000 

1111111111111111100000000 

1000000000000000000000000 
1111111111111111111111111 

1000000000000000000000000 
1111111111111111111111110 
1111111111111111111111111 

1000000000000000000000000 
1111111111111111111111111 

1000000000000000000000000 
1111111111111111111111111 

1000000000000000000000000 

1111111111111100000000000 

1100000000000000000000000 

1000000000000000000000000 
1111111111111111111111111 

1000000000000000000000000 

1111111111111100000000000 

1000000000000000000000000 

1111111111111111111111110 

1111111111111111100000000 
1111111111111111111111110 
1111111111111111111111110 
1111111111111111111111110 
1111000000000001111111100 
1111111111111111111111110 
1111111111111111111111111 


Scenario 20: WHTC with security. 


Hydro1 
Hydro2 
Hydro3 
Hydro4 
Hydro5 
Hydro6 
Hydro7 
Hydro8 


Daily operating cost:895407.39 $ 
Hours (0-24) 


1000000000000000000000000 
1000000000000000000011111 
1111111111111111111111111 
1000000000000000000000000 
1000000000000000000001111 
1000000000000000000000000 
1111111111111111111111111 
1000000000000000000000000 
1111111111111111111111111 
1000000000000000000000000 
1111111111111111111111111 
1000000000000000000000000 
1111111111111111111111111 
1000000000000000000000000 
1000000000000000000001111 
1000000000000000000000000 
1111000000000000000000000 
1111111111111111111111111 
1000000000000000000000000 
1100000000000000000000000 
1111000000000000000000000 
1111111111111111111111111 
1000000000000000000000000 
1111111111111111111111111 
1000000000000000000000000 
1111001111111111111111110 
1111111111111111111111110 
1111001111111111111111110 
1111111111111111111111010 
1111111111111111111111111 
1000000010000000011111110 
1000001101110111100101011 
1111111111111111111111111 
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Tables 2 and 3, in which 1 and O represent ON/OFF states of units 
at different hours, and hour 0 represents the initial condition. 

According to simulation of scenario 10 and 20, the following 
results can be obtained: 

In scenario 10, the thermal units 2 and 3 which have the high 
production cost is remained off on all period planning. While in 
scenario 20, both these units in the final of the day hours (from 
hour 20 to 24) have been forced to turn on and connected to 
the grid. 

In scenario 10, the thermal plant 40 is not committed at hour 
2-24, while in scenario 20 this plant is committed. In scenario 10, 
the thermal plant 44 is ON at total of 24 h while in scenario 20 it is 
ON only at hour 1-3. The expensive units of plant 52 is committed 
at hour 1-13 in scenario 10 while it will be ON at total of 24h in 
scenario 20. Also the expensive plant 7 and 39 has been ON more 
hours in scenario 10 than scenario 20. So expensive plant 7 is 
connected to grid from hour 1 to 16 while in scenario 20 it has 
been committed only at hour 21-24. 

The hydro units 6 and 7 in scenario 10 are committed more 
hours than scenario 20. 

The whole generating power of hydro and thermal plants in the 
scenarios 10 and 20 are 15107.94/54218.871 MW and 14656.002/ 
62380.987 MW respectively. As can be seen, in scenario with 
lower production cost, i.e. scenario 10 the amount of hydro plants 
production is more. Also the ratio of the hydro generating power 
to total of power generated with hydro and thermal plants is 
21.79% and 19.02% in scenario 10 and 20 respectively. 

In scenario 10, the thermal plants generate the 54218.871 MW 
power in all period of time and the average of generating power 
with those is 2259.119 MW. While the maximum generating 
power of thermal power is 3369.614 MW in hour 15. It is 
noteworthy that variation of thermal plants generating power is 
less compared with variation of hydro plants generating power 
given in next paragraph ((3369.614—2259.119)/2259.119=49.16%). 
The reason of these low variations is the limitation of thermal 
plants power variations and due to this constraint, the thermal 
units do not follow the rapid load variations. 

In scenario 10, the hydro units in all period of time produce the 
15107.94 MW that its average is the 629.497 MW at each hour. 
The minimum generating power at hour 14 is 104.415 MW and 
the maximum of hydro units generating power at hour 15 is 
1393.809 MW. The variations of hydro units generating power are 
very high and pursue the load variations in instantaneous market. 
Also the amount of whole SO and NOX gases produced during 
24 h has been 80697.95 and 41878.83 in scenario 10 and 101391.11 
and 50755.54 in scenario 20 respectively. This indicates that 
despite production cost in scenario 10 is less than scenario 20, 
the amount of pollution gases generation is less too. Also it should 
be noted that the minimum pollution gases (SO and NOX) 
generation related to scenario 10. 


Case 2. Wind-hydrothermal coordination without considering 
of security constraints 


Again each 20 scenarios have been simulated without con- 
sidering security constraints and their power production cost has 
been compared with each other. It can be seen that with not 
considering of the security constraints the minimum and max- 
imum of production cost daily related to the scenario 1 and 20 
with 698609.36 $ and 811174.28 $ respectively. The commitment 
schedule of scenario 1 is shown in Table 4. Also it is obvious which 
considering of security constraints will be caused to increasing the 
power production cost daily compared to when they are neglected. 
The main reason is that if prior due to satisfying security 
constraints some cheap units cannot be ON, now in this state 
can be ON. With comparing of scenario 1 (scenario with lowest 


Table 4 
Scenario 1: WHTC without security. 


Daily operating cost: 698636.98$ 
Hours (0-24) 


1-4 1000000000000000000000000 
5 1111111111111111111111111 
6-10 1000000000000000000000000 
11 1111111111111111111111111 
12-19 1000000000000000000000000 
20-21 1111111111111111111111111 
22-23 1000000000000000000000000 


1111111111111111111111111 
26 1000000000000000000000000 


27-29 1111111111111111111111111 
30-38 1000000000000000000000000 
39-40 1111111111111111111111111 
41-44 1000000000000000000000000 


45 1111111111111111111111111 


46-54 1000000000000000000000000 
Hydro1 1111111111111111111111110 
Hydro2 1111111111111111111111110 
Hydro3 1111111111111111111111110 
Hydro4 1111111111111111111111110 
Hydro5 1111111111111111111111111 
Hydro6 1000000000000000111110000 
Hydro7 1111111111111111111110000 
Hydro8 1111111111111111111111111 
Table 5 


Congested/uncongested status of violated lines in 
scenario 1 (with minimum cost) of Case 2. 


Hours (0-24) 


L37 0000000000000000011010000 
L41 0000000000000000000110000 
L124 0000000010001001011100000 
L126 0000000000001001000100000 
L127 0000000000001000000000000 
L134 0000000010001001011110000 
L136 0000000010001001011110000 
L137 0000000000001001000100000 
L138 0000000000001000000100000 
L139 0000000000001000000100000 
L142 0000000000000000111110000 


cost related to state without considering of the security con- 
straints) and scenario 10 (scenario with lowest cost related to 
state with considering of the security constraints), it is seen that 
the both inexpensive units 4 and 5 are committed in scenario 10, 
while in scenario 1 only unit 5 maintains its ON state. Also it is 
obtained that the both expensive thermal units 7 and 52 in 
scenario 10 to supplying load and satisfying of security constraints 
has been forced to be ON at hour 1-16 and hour 1-13 respectively, 
whereas in scenario 1 these unit all 24 h has not been committed. 
Compared scenario 1 with scenario 10, thermal unit 40 that was 
OFF in all of 24h has been ON and unit 39 is committed at the 
more hours to compensate for the reduced supply for decommit- 
ting of thermal units 19, 44 and 52 and to satisfy physical 
constraints in scenario 1. 

With comparing of scenario 20 (scenario with highest cost) 
with and without considering of the security constraints, it is seen 
that the economical unit 4 and both inexpensive units 19 and 52 
had been forced to be ON over a day to supplying load and 
satisfying of security constraints, have been OFF. Instead of these 
decommitted units, inexpensive unit 10 that was OFF before will 
be ON at total of 24 h. Also expensive thermal units 2, 3 and 7 and 
units 34 and 44 that were ON only in the first or last hours of the 
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Table 6 
Congested/uncongested status of violated lines in 
scenario 20 (with maximum cost) of Case 2. 


Hours (0-24) 


L30 0000000000000000000111000 
L37 0000000100000000000101010 
L41 1000001111000000001101111 
L54 0000000000000000000001000 
L124 0000001101000000000001011 
L126 0000000001000000000000000 
L127 0000000000000000000000000 
L134 0000001101000000000101010 
L136 0000001101000000000101010 
L137 0000000000000000000000000 
L138 0000000000000000000000000 
L139 0000000000000000000000000 
L142 0000000000000000001101110 


Table 7 
Scenarios with minimum/maximum daily operating cost ($) in 3 cases. 


With security constraints Without security constraints 


Case 1 Case2 Case 3 Case 1 Case 2 Case 3 
Scenario 10 731831 -----  ------ 9 ------ ------ 9 ------ 
Scenario 20 89540738 ---- ——=—-— 811174.28 ------ 
Scenario 1 ------—  ------ ss = 698609.36 ------ 


Scenario 15 | ---- 


Table 8 
Aggregated results of the stochastic security-constrained WHTC framework in 
3 cases. 


Case 1 Case 2 Case 3 


Daily operating cost ($) 788596.71 730831.36 1086309.34 


day, and also without considering security constraints during 
these hours they have been OFF and decommitted over a day. 
For compensating of reduced supply due to decommitting these 
5 units, unit 39 that was ON only in first 3 h of day it have been 
forced to be ON over a day to supplying load. 

When in this case, we calculate the WHTC solution by excluding 
transmission and voltage constraints, transmission flow violations in 
scenario 1 and scenario 20 occur in different lines and hours. 
Tables 5 and 6 show all of transmission flow violations on congested 
lines for these scenarios without considering security constraints, in 
which 1 and O represent congested/uncongested status of lines at 
different hours, and hour 0 represents the initial condition. 

Results show that branch 41 is very congested. This branch has 
80 MW of capacity and that is not sufficient to transmit less 
expensive generation from the right-hand side of the system to 
the left-hand side and the its’ transmission flow in Case 3 is near 
or at the capacity limit at hours 1, 2, 8, 9, 10, 11, 12, 14, 15, 16, 18 and 
19. Also scenario 20 has less congested line than scenario 1, for 
example congested lines 127, 137, 138,139 in scenario 1 do not have 
any violated flow in scenario 20. 


Case 3. Evaluation of effect of hydro units on 
the security-constrained WHTC problem 


Once again, each 20 scenarios have been simulated and 
investigated. But this time with eliminating of hydro plants from 
problem and then their power production cost that related to the 
thermal and wind units, are compared with each other. In this case 


the minimum and maximum production cost daily are related to 
the scenario 15 with 991005.67 $ and scenario 20 with 1204288.30 
$ respectively. Because of not produce power by hydro units which 
are inexpensive units, it can be seen that production cost in each 
scenario is higher than two previous cases included hydro units. 
For quick reference, the three case tests are briefly described in 
Table 7. Also, the aggregated results of the stochastic security- 
constrained wind hydrothermal coordination framework, accord- 
ing to the scenario aggregation procedure, are given in Table 8. 
Using the proposed stochastic framework, all 20 accepted scenar- 
ios contribute into determining the security-constrained wind 
hydrothermal coordination results according to their probability 
values. 


5. Conclusion 


In this paper, the thorough and comprehensive MIP formula- 
tions for solving security-constrained wind hydrothermal coordi- 
nation problem at day-ahead are proposed. The objective is to 
minimize the daily operating cost of ISO. The proposed framework 
consists of many practical constraints such as prohibited operating 
zones, valve loading effects, crew constraints, dynamic ramp rate, 
and fuel and emission limitations for thermal units. Furthermore, 
for the hydro plants, multi performance curve with spillage and 
time delay between reservoirs are considered. We also consider 
transmission network constraints. The main feature of the pro- 
posed framework refers to the linear nature of the formulations 
which is very important for application of the model in the large 
scale and the real size power system. Therefore, the proposed 
scheme is practical to generate appropriate information for ISOs. 
Roulette wheel mechanism and Monte Carlo Simulation (MCS) are 
employed for random scenario generation wherein the stochastic 
optimization problem is converted into its respective deterministic 
equivalents. Results of the implementing the proposed model 
show that the solution time of the problem is rationale. The 
research work under way in order to reduce the computational 
burden due to increasing the problem size, two approaches (and 
their combinations) can be taken. These are the use of decom- 
position approaches such as benders decomposition and parallel 
computing used for each scenario. 
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Appendix A 


Auto and cross correlation coefficients related to two imple- 
mented wind farms [10]: 


w= O11 O12 Se a) 
1 \ ða 02] \-044 0.22 


jz o1 \ _ / 0.008 
1 \ 002 a naar) 
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